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ABSTRACT 

Analysis of the Cassini Ultraviolet Imaging Spectrometer (UVIS) stellar and 
solar occultations at Titan to date include 12 species: N 2 (nitrogen), CH 4 
(methane), C 2 H 2 (acetylene), C 2 H 4 (ethylene), C 2 H 6 (ethane), C 4 H 2 (diacety- 
lene), CeH 6 (benzene), CeN 2 (dicyanodiacetylene), C 2 N 2 (cyanogen), HCN (hy- 
drogen cyanide), HC3N (cyanoacetylene), and aerosols distinguished by a struc- 
tureless continuum extinction (absorption plus scattering) of photons in the EUV. 
The introduction of aerosol particles, retaining the same refractive index proper- 
ties as tholin with radius ~125 A and using Mie theory, provides a satisfactory 
fit to the spectra. The derived vertical profile of aerosol density shows distinct 
structure, implying a reactive generation process reaching altitudes more than 
1000 km above the surface. A photochemical model presented here provides a 
reference basis for examining the chemical and physical processes leading to the 
distinctive atmospheric opacity at Titan. We find that dicyanodiacetylene is con- 
densable at ~650 km, where the atmospheric temperature minimum is located. 
This species is the simplest molecule identified to be condensable. Observations 
are needed to confirm the existence and production rates of dicyanodiacetylene. 

Subject headings: planetary systems — radiative transfer — atmospheric effects- 
planets and satellites: individual (Titan) — methods: data analysis, numerical 

1. Introduction 

Titan is Nature's laboratory for organic synthesis. The major molecules in the atmo- 
sphere are N 2 and CH4. The coupled chemistry between nitrogen and carbon leads to high 
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abundances of nitrogen/carbon compounds, such as hydrogen cyanide (see, e.g.. lYung et al. 



1984 ; ICoustenis et al.l 120031 ; IWilson fc Atreyal 120041 ) . This is c aused primarily by the low 



gravity of Titan which allows hydrogen to escape readily (e.g., lYelle et al.l 120061 ) . resulting 
in low hydrogen abundance and rich hydrocarbon production. When the order of hydro- 
carbons and nitriles is large enough, they condense to form aer osols and precipitat e. The 
stratosphere is a major region for the production of haze (e.g., McKay et al.l l200ll ) . The 
estimated vertically integrated rate is 0.5-2 xl0~ 14 g cm -2 s~ x with h aze formation taking 
place in the 300 km to 500 km region. Total mass loading according to iMcKay et al.l (1200 ll ) 
was about 250 mg m -2 , made by molecules with a C/N ratio of 2-4 and a C/H ratio of about 
unity. The present work obtains about 100 mg m -2 (assuming density 3 g cm -3 ) above 300 
km. 



2. Cassini UVIS Observations 

On 13 December 2004, the Cassini UVIS recorded the occultation of two stars, A Sco 
(Shaula; latitude —36°) and a Vir (Spica; latitude ran ge of +63° and +48°), near the end 



of th e second Titan flyby (T r) (details are referred to IShemansky et al.l 12005c IShemansky 



20061 ; IShemansky et al.l 120071 ) . The fully reduced results from A Sco only are referenced 
here. The vertical distribution of the aerosol component as the terminal product of N 2 /CH 4 
physical and chemical processes, is of primary interest to this Letter. The spectral region 
1850-1900 A (SP1) is effectively free of the measurable hydrocarbon and cyano species in 
this atmosphere apart from aerosols, and this region is used photometrically to trace the 
aerosol structure. The vertical profile of extinction for SP1 below 1000 km was carried out 
at the highest possible ray-height resolution (3-5 km) in order to reveal possible structure in 
the distribution. 

The best fit to the transmission spectrum at impact parameter h = 514 km is shown 
in Figure [1] as an example of spectral reduction. Table [1] shows the extracted line-of-sight 
abundances from this element of the occultation. Figure [2] shows the aerosol density vertical 
distribution (heavy dots) derived from the reduction of the A Sco occultation, compared to 
the derived CH4 profile (dashed line). The reduced data extend from h = 330 km to 970 km 
where signal noise terminates the reduction (see Shemansky et al. 2007 for details). The 
remarkable property of the aerosol distribution is the sudden departure from tracking the CH 4 
abundance at h = 468 km toward higher altitudes. The interval between 468 km and 550 km 
where abundance remains approximately constant, is interpreted as a primary aerosol source 
region. The density extraction is a direct deconvolution of the abundance distribution. The 
derivation is based on assumed spherical uniformity and assumed uniformity in composition 
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with altitude. The observations show measurable simple scattering in the 1800 A region at 
h = 1250 km, indicating that aerosols are extensively distributed into the thermosphere. 

The cross sections for extincti on of UV light by a erosols were computed using the com- 
plex refractive index measured by iKhare et al.l (119841 ) for solid state tholins, and the scat- 



tering code of Mishchenko fc Travisf 



19981 ). Representative complex indices of refraction at 



588, 1215, 1631 and 2384 A are (0.963, 0.62), (1.74, 0.37), (1.65, 0.24) and (1.68, 0.21), 
respectively. Assuming a mean radius of 125 A for aerosols, the corresponding extinction 
cross sections are 2.4xl0~ 12 , 3.9xl0~ 13 , 2.0xl0~ 13 and 1.2xl0~ 13 cm 2 , respectively. The 
actual computation used in this work has a finer wavelength grid. Note that the current 
UVIS data set no constraint on the vertical variation of the optical property of aerosols. 

Extinction of the EUV stellar photons is dominated by CH4 in the 1100 A to 1400 A 
region. At longer wavelengths the structure is a combination of the higher order hydrocarbon 
and cyano species. The C 6 H 6 cross section peaks at 1759.9-1815.1 A, blended primarily with 
C 2 H 4 . Dicyanodiacetylene has a cross section peak in the SP1 spectral r egion. Dicyanodi- 
acety lene and benzene have not been detected in the absorption spectra (iShemansky et al. 
20071 ). Aerosol extinction is detectable at ~970 km in the transmission spectra and domi- 
nates all absorbers at all wavelengths in the UVIS except for CH 4 , at altitudes below 400-450 
km. 

Figure H] shows the contribution of aerosol extinction to the total measured optical depth 
at 514 km. In the spectral region SP1 aerosol extinction is entirely responsible for the optical 
depth. The measurable spectral region for extracting the aerosol component is 1500-1900 
A, where the wavelength dependence of extinction shows a proportion ality to A~ 1,5 . The 



Voyager and Cas sini photometric observations in the UV spectral region (IPorco et al. 112005 



West et al.ll2006l ) have revealed the presence of detached haze layers at Titan. The Cassini 
results show the presence of a latitudinally uniform detached layer near 500 km in forward 
scattered 3380 A photons. The relationship of this phenomenon to t he aerosols identifi ed here 
requires further investigation. A comparison to Voyager results ( jSmith et al.l [1982J ) shows 
that the major differences are that the apparent strong extinction by aerosols takes effect 
about 100 km higher for Voyager and significantly more extinction is evident for Voyager in 
the 700-1000 km region. The Voyager data show a broad extinction maximum near 770 km. 



3. Photochemical Modeling 



Vertical profiles of the major species h ave been calculated using a phot oche mical model. 
The photochemical reactions are taken from lYung et al.l (119841 ) . lYungi (119871 ). and lMoses et al 
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( 12000I ). The chemica l scheme to C6N2 is hypothesized to be similar to that to C4N2, as de- 
rived by mind (119871 ) and summarized in Table [2j The temperature profile is based on the 
Cassini measurem ents (Figured], dotted line). The vertical eddy mixing profile is taken from 
Yung et al.l (Il984f ). The model simulation is diurnally averaged at low latitude. The incident 
UV flux is the mean between solar maximum and minimum. Table [TJ provides a summary of 
model results. Sensitivity to the selection of hydrocarbon kinetics and that of kinetics and 
vertical eddy coefficients are shown by models D and WA04, respectively. 

We fix the N2 abundance to that derived from the Cassini measurements. The model 
starts with a hydrostatic atmosphere. With the prescribed vertical diffusion coefficients and 
taking the photolysis of CH 4 into account, the abundance of CH 4 is overestimated (see Ta- 
ble [TJ), compared with the measurements. To bring the model into better agreement with 
the observations, we introduce an ad hoc advection which transports species other than N2 
downward. The wind is prescribed with strength proportional to the inverse of the square 
root of atmospheric density. The wind speed reaches -20 cm s -1 at the top of the model at- 
mosphere (~1500 km). The assumed downward wind is qualitatively consi stent with global 



circu lation that has a downward transport at mid to high latitudes (e.g., iLebonnois et al. 
20011 ) . A comparison with modeled CH 4 abundances between models A and C is shown 
in Figure [3j We note that dynamics plays an important role in distributing photochemical 
products (e.g.. ILebonnois et al.ll200ll ). especially in the regions above ~500 km (the regions 



of interest to this work) where the transport time is, in general, shorter tha n the chemical 
removal time of hydrocarbon and cyano species (e.g., 



Lebonnois et al. 


2001). 


(e.g., 


Flasar et ali 


2005) 



Wilson &: Atreyall2004l ); current simu- 
imited to the regions below ~400 km 



20051 ) are the consequence of atmospheric transport and photochemical 



processes. 

The modeled profiles of HCN, HC 3 N, C 6 H 6 , and C 6 N 2 are presented in Figure [3J Cal- 
culated abundances are compared to extraction from observation at 514 km in Table [TJ 
Comparisons between models and observations at other impact parameters will be deferred 
to a later paper. The results for five variations on the model at this impact parameter are 
given in Table [TJ In general, our base models (models A and B) overestimate the abundances 
of hydrocarbons by as much as 10 times. An indication of the difference in the predicted 
model C and observed optical depth spectra is shown in Figure [TJ Model C (Table [TJ is 
too high relative to measured abundance in C2H6 and C4H2. The modeled C6N2, CeHg, and 
HC3N are also well above the upper limits set by observation. There are two ways of reducing 
the abundance s. (1) Faster transport of photochemical products to the lower atmosphere as 
in the model of I Wilson Atreyal (120041 ) (model WA04). Comparing the transport time con- 
stant with the chemical destruction time, the abundances of C2H2, C2H6, HCN, and C2N2 



- 5 - 



have sensitivity to transport and those of C2H4, C4H 2, CfiN2, CfjHfj, and HC 3N at ~500 km 



are close to being in photochemical equilibrium (e.g.. IWilson fc Atreyal 120041 ). (2) Relatively 
rapid two-body physicochemical processes forming aerosols as simulated in models C and D. 
The loss rates in these models are assumed to be proportional to the physical collision rates 
between aerosols (with radius 125 A) and molecules; adsorption reactions are assumed for all 
photochemical species listed in Table [TJ The aerosol density is from Figure [2j The adsorption 
efficiency for this absolute loss, an assumed value of 0.01, yields the concentrations shown in 
Table [TJ As described below, this process is required for explaining the aerosol abundance 
shown in Figure [2j 



4. Discussion 

The source of aerosols has long been a puzzle in the atmosphere of Titan. It is generally 
believed that the synthesis of increasingly complex hydrocarbon and nitrogen compounds 
will eventually lead to saturation, resulting in coagulation and precipitation. However, the 
chemical composition of the condensible species has not yet been established. In this Let- 
ter, we propose that the simplest condensible compound is C6N2. The abundances of the 
higher order species in the UVIS observations are significantly lower than the present model 
calculations. Assuming that the model conversion rates for the CH4 source are basically 
correct, there is an implied loss rate for these species that is substantially higher than the 
model provides. The model calculation contains an absolute loss to the measured aerosols 
using a conservatively small adsorption probability. If a significant fraction of the implied 
loss is delivered to the production of aerosols, the model can be adjusted by assuming a 
larger irreversible adsorption probability and a consequent higher precipitation rate for the 
aerosols. This will bring the model abundances into better conformance with observation, 
but will not necessarily resolve differences in partitioning, and may not resolve issues raised 
in regard to rate limits for the mass flow process that cycles to the surface. We propose that 
the aerosol formation is initiated by the condensation of C6N2 and adsorption to external 
meteoritic dust; both serve seeds for aerosol formation. The subsequent physical processes 
of adsorption on the existing aerosols and (photo) chemistry converting these clusters into 
refracto r y tho lins constitute the production and maintenance of the aerosol distribution. 



Huntenl (120061 ) has recently proposed that haze is a major sink of ethane at Titan. The pro- 
cess of formation of aerosols requires the stable adsorption of the higher order hydrocarbon 
and nitrile species. 



The production rate of C6N2 in the model limits the production rate of aerosols from 
this direct path. The saturation density of dicyanodiacetylene shown in Figure [3] shows 
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that condensation can take place between ~550 and 800 km. The volume production rate 
of CeN 2 in this region is quite uniform (~1 molecules cm -3 s -1 ); the column integrated 
(550-800 km) rate is ~10 7 molecules cm -2 s _1 , or ~2xl0~ 15 g cm -2 s _1 . This contribution 
is a small fraction of the total aerosol production, but is extremely important as a source 
of condensation nuclei. The maximum rate of aerosol production, however, is set by the 
total photolysis rate of CH4 which is on the order of 10 10 molecules cm" 2 s _1 , or 2xl0~ 13 
g cm" 2 s" 1 of ca rbon, a result that is independently corroborated by the H 2 escape flux 



(lYelle et al.ll2006l ). The rate of removal of molecules by the existing particles in model C is 
close to this maximum rate. However, the resulting aerosol profile (thin solid line in Figure 
2) underestimates the observed aerosol abundances above ~800 km, suggesting an additional 
source at the top of the atmosphere, equivalent to a downward flux of 5xl0~ 14 g cm -2 s" 1 
(thick solid line). This flux is consistent with th at inferred from the Cassini Ion Neutral Mass 



Spectrometer measurements (IWaite et al.ll2007l ) on the basis of ion chemistry not considered 



here. Note that a sedimentation velocity of 0.25 cm s -1 has to be imposed in order to match 
the observed aerosol profile; this implies that the radius of the aerosols must be ~125 A in 
the mesosphere and thermosphere. 

We emphasize that there are significant uncertainties in rate process quantities and 
stronger constraints are needed. Laboratory measurements for the adsorption rates on 
aerosols in collision with the high order molecules are required to verify this process, and 
to provide a better constraint to the aerosol mass loading in the atmosphere of Titan. In 
addition, the photochemical paths to high order hydrocarbon and nitrile compounds such as 
CgN 2 and CqHq are speculative. Atmospheric dynamics in the mesosphere of Titan also plays 
a central role in molecule/aerosol mixing. Further laboratory measurements and Cassini ob- 
servations will provide valuable information in refining our understanding of the chemical, 
dynamical, and microphysical processes in the atmosphere of Titan. 
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Fig. 1. — The transmission spectrum of the UVIS A Sco occultation integrated over 
the impact parameter 514 km to 537 km (light line), compared to the best-fit synthesis 
(jShemansky et al.ll2007l ) using the combined identified species (heavy line), and to model C 
from the present physical chemistry code (dashed line). The dotted line shows the aerosol 
component combined with the CH 4 absorber in this reduction. The CH 4 absorber (A < 
1490 A) is included with aerosol here as a means of including a large part of the impact of 
instrument point spread function on the fitting process. The optical depth in region SP1 
is entirely attributed to aerosol extinction; The small difference between the observed data 
and the aerosol component at SP1 is an artifact of the UVIS EUV instrument point spread 
function. The abundances of the species for this case are given in Tabled] See text. 
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Fig. 2. — Aerosol density (heavy dots) der ived from the UVIS A Sco occupation com pared 
to the CH 4 (dashed line) scaled by 10~ 9 (jShemansky et al.ll2005l : IShemanskyl 120061 ). The 
increase of the mixing ratio of the UVIS aerosols through the mesosphere to at least 1000 
km implies that the production of aerosols must take place at significant rates throughout 
the mesosphere and thermosphere. The UVIS derived temperatu re profile is shown by th e 



dotted line, which reflects a correction to the one presented by IShemansky et al.l (120051 ) . 



The model aerosol profiles are shown by the thin and solid lines (see text). The over- and 
underestimations are due to the fa ct that we assume a constant sedimentation velocity of 
0.25 cm s _1 , calculated at 535 km (ICabane et al.lll992l ). 
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Fig. 3.— Modeled (model C) vertical profiles for CH 4 (thick solid), HC 3 N (dashed), HCN 
(dash-dotted), C 6 N 2 (triple-dot-dashed), and C 6 H 6 (long-dashed). Thin solid line represents 
modeled CH 4 by model A. The s aturation densit y of C 6 N 2 extrapolated from high temper- 
ature measurements (295-369 K. ISaggiomolll957l ) is shown by dotted line. The resulting H 2 
(3xl0~ 3 ) and CH4 (2.3%) mixing ratios at 1174 km and H 2 escape flux (7xl0 9 molecules 
cm -2 s" 1 ) at the top are in good agreement with th e observations (4 ±lxl0~ 3 , 2.7±0.1%, 
and 1.2±0.2xl0 10 molecules cm 2 s , respectively) (lYelle et al.ll2006l ). 



Table 1. Summary of Model Results 



Molecule 


Cassini 


ivioaei a 


ivioaei 


Model u 


Moflei L) 


WAU4 


N 2 (xlO 21 ) 


5.8 


5.8 


5.8 


5.8 


5.8 


5.8 


CH 4 (xlO 19 ) 


6.0 


15 


9.4 


9.5 


9.7 


13 


C 2 H 2 (xlO 17 ) 


2.1 


15 


9.1 


1.9 


1.7 


1.5 


C 2 H 4 (xlO 16 ) 


4.0 


9.3 


5.7 


4.0 


2.0 


3.4 


C 2 H 6 (xlO 16 ) 


7.0 


200 


110 


17 


9.2 


20 


HCN (xlO 17 ) 


1.0 


5.6 


3.7 


0.69 


0.53 


0.017 


C 4 H 2 (xlO 15 ) 


4.5 


59 


37 


12 


2.1 


41 


C 6 N 2 (xlO 14 ) 


<1.0 


15 


16 


5.7 


8.0 




C 6 H 6 (xlO 14 ) 


<1.4 


18 


13 


12 


0.041 


0.21 


HC 3 N (xlO 15 ) 


<3.9 


27 


25 


6.9 


7.8 


0.96 


C 2 N 2 (xlO 15 ) 


<4.0 


0.84 


1.2 


0.34 


0.42 


0.00035 


Tholin (xlO 11 ) 


4.6 













Note. - Values are line-of-sight column integrated abundances, in molecules 
cm -2 , reported by matching the observed N 2 abundance. Model A: hydrostatic 
atmosphere. Model B: non-hydrostatic atmosphere, an ad hoc downward wind 
and extinction due to the derived tholins are assumed (see text). Model C: same 
as Model B but also with additional sinks for the tabulated nine photochemical 
species (see tex t). Model D: same as model C but with the upd ated hydrocarbon 
chemi stry from lMoses et al.l (120051 ) . WA04: model results from lWilson & Atreya 
(120041 ) . Note that in this Letter, the microphysical processes of CgN 2 are not 
solved self-consistently; and hence, the tabulated abundances of CeN 2 do not 
reflect the removal by condensation. 
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Table 2. Chemical Reactions to C6N2 



Label 


Reactants 


Products 


Rate Coefficients 3 


R454 


HC 5 N + mv 


-> C 4 H + CN 


= J(HC 3 N + hu - 


4 C 2 H + CN); (1) 


R455 


HC 5 N + mv 


-> C 5 N + H 


= J(HC 3 N + hu - 


4 C 3 N + H); (1) 


R456 


C 6 N 2 + hu 


-> C 5 N + CN 


= J(C 4 N 2 + hu -> 


C 3 N + CN); (2) 


R492 


C 3 N + HC 3 N - 


-> C 6 N 2 + H 


= k(C 2 H + C 2 H 2 


-> C 4 H 2 + H); (3) 


R495 


CN + C 4 H 2 


-> HC 5 N + H 


= k(CN + C 2 H 2 - 


-> HC 3 N + H); (2) 


R496 


CN + HC 5 N - 


-> C 6 N 2 + H 


= k(CN + HC 3 N - 


^ C 4 N 2 + H); (2) 


R497 


C 5 N + CH 4 


-> HC 5 N + CH 3 


= k(C 3 N + CH 4 - 


-> HC 3 N + CH 3 ); (2) 


R498 


C 5 N + C 2 H 6 - 


-> HC 5 N + C 2 H 5 


= k(C 3 N + C 2 H 6 


^ HC 3 N + C 2 H 5 ); (2) 



a Estimated from the quoted reactions. Units are s -1 for photolysis reactions (J) and 
cm 3 s _1 for two-body reactions (k). The p hotolysis rate coefficien ts are given at the top 
of the mo del atmosphe re. References: (1) 



((200J); (2) lYungl (119871 ): (2) iLebonnois et al. 



Lebonnois et al. and IWilson fc Atreval 



(120011) and lOpanskv fc Leonel jl99fih 



